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We consider diffusive iattice gases on a ring and anaiyze the stabiiity of their density profiles 
conditionally to a current deviation. Depending on the current, one observes a phase transition 
between a regime where the density remains constant and another regime where the density becomes 
time dependent. Numerical data confirm this phase transition. This time dependent profile persists 
in the large drift limit and allows one to understand on physical grounds the results obtained earlier 
for the totally asymmetric exclusion process on a ring. 
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I. INTRODUCTION 

Stochastic lattice gas models have been exten- 
sively studied recently as they are among the sim- 
plest examples of non-equilibrium systems. A pow- 
erful approach to understand their steady state 
was developed by Bertini, De Sole, Gabrielli, Jona- 
Lasinio, Landim, as a macroscopic fluctuation the- 
ory (MFT) which gives, for large diffusive systems, 
the probability distribution of trajectories in the 
space of density profiles 0, Q . The MFT relies on 
the hydrodynamic large deviation theory 0- [J 
which provides estimates for the probability of ob- 
serving atypical space/ time density profiles. It gives 
a framework to calculate a large number of proper- 
ties of stochastic lattice gas, such as the large devia- 
tion functional of the density profiles. Recent devel- 
opments of the hydrodynamic large deviation theory 
0, Q enabled to estimate also the large deviations 
of the current through the system. What the MFT 
provide are the equations of the time evolution of 
the most likely density profile responsible of a given 
fluctuation. What it does not provide, in general, 
is the solution of these equations which would give 
quantitative predictions for the distribution of the 
fluctuations. So far, for the large deviation func- 
tional of the density profiles in the steady state, the 
equations could only be solved in a few cases of non- 
equilibrium systems with open boundaries (the sym- 
metric exclusion process 0, the Kipnis, Marchioro, 
Presutti model [1 ®). For the SSEP the results 
of the MFT were in full agreement with the results 
obtained [ToL ITTl IT^ from the exact knowledge of 
the weights of the microscopic configurations in the 
steady state. 

In our previous work , we developed a theory to 



calculate the large deviation function of the current 
through a long one dimensional diffusive lattice gas 
in contact at its two ends with two reservoirs at un- 
equal densities. Our approach was based on an as- 
sumption, the additivity principle, which relates the 
large deviation function (LDF) of the current of a 
system to the LDF's of subsystems, when one breaks 
a large system into large subsystems. This assump- 
tion is in fact equivalent to the hypothesis, within 
the hydrodynamic large deviation framework, that 
to observe, for a very long time period T, an aver- 
age current q = Qt/T, the system adopts a profile 
with a shape, fixed in time, but of course depend- 
ing on q (here Qt is the total number of particles 
transfered, say from the left reservoir to the system 
during time T). The additivity principle allows one 
to obtain explicit expressions \J\ for all the cumu- 
lants of the integrated current Qt- The predictions 
of our theory were tested in a few cases pfl and the 
results were found in complete agreement with what 
was already known or what could be derived by al- 
ternative approaches [l3L IliL Il5| . 

Recently, it was pointed out [|| that even if our 
predictions are valid for some diffusive lattice gas, 
it might happen that, to produce an average current 
q over a long period of time, the best profile is time- 
dependent. One of the goals of the present work 
is to show that, for a simple example, the weakly 
asymmetric exclusion process on a ring, this indeed 
happens for some range of parameters. 

Let us consider, as we shall do it in the rest of 
this paper, the time evolution of a one dimensional 
stochastic lattice gas on a lattice of N sites. Accord- 
ing to the hydrodynamic formalism a given lat- 
tice gas can be characterized by two functions D{p) 
and <j{p) of its density p. One way to define them 
is to consider a one dimensional system of length 
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N connected to reservoirs at its two ends. For such 
a lattice gas, the variance of the total charge Qt 
transfered during a long time T from one reservoir 
to the other is given, for large N, by definition of 
a(p) by 



(Qt) 
T 



N 



fl 



when both reservoirs are at the same density p. On 
the other hand if the left reservoir is at density p + 
Ap and the right reservoir at density p, the average 
current is given, for small Ap, by 



(Qi 



T 



N 



(2) 



which is simply Fick's law and defines the function 
D(p). In the symmetric simple exclusion process 
a(p) = p(l — p) and D(p) = 1/2 whereas in 
the Kipnis, Marchioro, Presutti model <r(p) = 
p 2 and D(p) = 1/2. The effect of a uniform weak 
electric field of strength v/{2N) acting from left to 
right on the particles is to modify (J2J into 



(Qt) D(p)Ap , va{p) 



T 



N 



N 



(3) 



This equation follows from the linear response the- 
ory 0. 

Once D(p) and cr(p) are known for a given diffu- 
sive system, the probability of observing the evolu- 
tion of a density profile p(x, s) and a rcscalcd current 
j(x, s) for < s < T during a time T ~ TV 2 is given, 
according to the hydrodynamic large deviation the- 
ory @, by 



Pio(j(x,s),p(x,s)) 



exp 



~^[0,T] w> p) 
N 



(4) 



where TL T j is denned by 

■£[o,t] (i> p) 



T ,1 

ds / dx 
o Ja 

[j(x, s) + D(p{x, s))p'(x, s) - va{p{x, s))f 



(5) 



2a(p(x,s)) 

with p' = dp/dx and where the rescaled current 
j(x, s) is related to the density profile p(x, s) by the 
conservation law 



dp(x,s) dj(x,s) 



ds 



dx 



(G) 



A formalism equivalent to this hydrodynamic large 
deviation theory was developed independently |l6l 
0] in the context of the full counting statistics of 
the transport of free fermions through disordered 



wires. A simple derivation of (JSJ) and © is given in 
the appendix. 

The large deviation function G(jo) of the current 
is then defined as 



Pro 



Or 
T 



3o_ 

N 



T 



for large T and N 



(7) 



(In (7J|, one has first to take the limit T — > oo and 
then make TV large; in practice J7J should hold when 
T ^ N 2 as N 2 is the characteristic time of a diffu- 
sive system of size N). 

Now according to Q), the large deviation function 
G(jo) is given by 



G(jo) 



lira 



1 



mil \ j [o,t] (i> p) 



(8) 



where the current j(x,s) satisfies for large T and all 
x the constraint 



i r 

1™ 7^ / j(x,s)ds =j 



(9) 



with the profile p(x, s) and the current j(x, s) con- 
nected by ©■ 

In the following, we will often consider, instead 
of Q), the generating function of the current: 

(e AQr ) - e T ^ for large T (10) 

and then, according to 17I8J1 . p(X) is given by 



p(X) = — max[Ajo + G(j )} 
iv jo 



(11) 



1 ,. 1 
— lim — max 

N T^oo T p(x.s) 



j(x,s)ds -I[ 0iT] (j,p) 



For a system of length N connected to two reser- 
voirs at densities p a and pb at its two ends, the 
calculation of the large deviation function G(jo) of 
the current is therefore reduced to finding the time- 
dependent profile p(x, s) which optimizes (JSJ under 
the constraints © and © and with the additional 
boundary conditions p(0, s) — p a and p(l,s) = pb- 
All the results of our previous work 7] follow then 
from the assumption that this optimal profile does 
not vary with time (except from boundary effects 
near time and time T which do not contribute in 
the large T limit). 

For a system on a ring of N sites, as we shall 
consider here, the optimization problem is the same 
except for the boundary condition which becomes 
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p(0,s) — p(l,s) and the fact that the total den- 
sity po on the ring becomes an additional conserved 
quantity 



p(x, s)dx = po 



Our paper is organized as follows: in section ITU 
we consider a general lattice gas on a ring and show 
under what conditions the flat profile becomes un- 
stable. In sections [TJ] and lllll we write the large 
deviation function of the current, when the optimal 
profile has a fixed shape moving at a constant ve- 
locity. In section IIVI we present exact numerical 
results on the weakly asymmetric exclusion process 
for small system sizes which give evidence that for 
some range of parameters, consistent with the re- 
sults of section ITTT1 the optimal profile is no longer 
flat but becomes space-time dependent (for p = 1/2 
it is only space dependent). In section we ana- 
lyze the limit of a strong asymmetry and obtain a 
simple expression for the large deviation function 
of the current, under the assumption that the opti- 
mal profile becomes in this limit a step function. In 
section fvTl we show that, even for a strongly asym- 
metric case such as the totally asymmetric exclusion 
process, one can exhibit time-dependent profiles de- 
termined by the Jensen Varadhan functional which 
give, for large system size, the exact large deviation 
function of the current previously calculated by the 
Bethe ansatz. 



II. CONDITIONS FOR THE STABILITY OF 
A FLAT PROFILE 

Consider a lattice gas on a ring of N sites with 
total density po. Under the assumption that the 
optimal profile p(x, s) is flat, that is 

p(x,s) = p a 

the large deviation function G(jo) is given 1~>IS[| by 

Gflat Jo) = j \ (12) 

2ct(/? ) 

which, by (JTTJ, gives for p(A) 

A(A + 2u)a(po) 



Mflat(A) 



2N 



(13) 



in which case due to © 

p(x,s) = po + — [-&(x)sai(ws) +j' 2 (x) cos(ws)] 
to 

where j\ (x) and j 2 {x) are periodic functions of pe- 
riod 1. The resulting expression for G(jo) to second 
order in e is 



G(jo) = - 



(jo ~ ^cr(Po)) 2 

2a(p ) 



f ■ i dx 
'o 



3i+32 
4ct(p ) 



(j'f+3^)D(Po) 2 , 3otii32-323[y(po) 
4w 2 cr(po) 2wcr 2 (p ) 

j 2 a"( Po ) u 2 a"(p ) j 2 a' 2 (p Q ) 



8w 2 (T(po) 2 



8cu 2 



4w 2 cr 3 (p ) 



As this expression is quadratic in the currents j\ 
and j2, the various Fourier modes are not coupled. 
Choosing for ji (x) and j 2 (x) 



ji(x) = acos(2irx) + &sin(27ra;) 
jiip) = ccos(27rx) + dsin(27ra;) 



one gets 



G(jo) 



(jo - va{po)f 



2a (po) 
-(a 2 + b 2 + c 2 + d 2 ) 



(ad — be) 
1 



(14) 



jo0-'(po)7T 



oj<j 2 (po) 
2ir i D 2 (p ) 



8a(p ) w 2 er(po) 
* 2 v' 2 (po)fo , v 2 -* 2 °"(Po) n 2 a"(p )j 2 ^ 



2^ 2 a3(p ) 



4lj 2 



Auj 2 a 2 (po) 



The flat profile is stable against the perturbation 
(|14(l . if this is a negative definite quadratic form in 
a, b, c, d. This is achieved when for all lu 
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2^D 2 (p ) | 7r 2 a' 2 (p )j 2 

8a(p ) uj 2 a(p Q ) 2u 2 <j 3 (p ) 



uW(p ) n 2 a"(p )j 2 



> 



joe'iPo)* 



2wcr 2 (po) 



4w 2 Alu 2 a 2 (p a ) 

The flat profile becomes therefore unstable if 
1 , 2^D 2 (p a ) , 7r 2 a' 2 (p )j 2 



8a(p ) LU 2 a(po) 2cj 2 cr 3 (p ) 

< 



v 2 -k 2 o"(pq) 7rV'(po)jo j cr'(po)TT 



Alo 2 (t 2 (po) 



2wer 2 (p ) 



(15) 



A natural question is whether one could increase 
G(jo) (and p(X)) by adding to this flat profile some 
small (e -c 1) space and time dependent perturba- 
tion of the form 

j(x, s) = j Q + e[ji(x) cos(ws) + j 2 (x) sin(a;s)] 



87T 2 D 2 (p )a(p ) + (v 2 cr 2 (po)-j 2 )a"(p Q ) < (16) 
which, given I|11I12I13|I . can be rewritten as 

47r 2 D 2 (p Q ) < iVp flat (A)a"(p ) (17) 
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When the flat profile becomes unstable, according 
to (|15|l . the current takes the form 

j(x, t) =j + A cos 2tt (x - x - J °J. ^ t\ (18) 

where the amplitude A would be determined by ex- 
panding G(jo) to higher order in e. 

One could analyze in a similar way the stabil- 
ity of the flat profile against other modes by choos- 
ing ji(x) = acos(27mx) + 6sin(27rnx) and j 2 {x) — 
ccos(27rna;) + ds'm(2irnx) and the threshold 1)16(1 
would become 

Sir 2 D 2 (p )n 2 a(p ) + (v 2 a 2 (p ) - jl)a"(p ) < 

This shows that that the fundamental (n = 1) is the 
first mode to become unstable. 



III. A SIMPLE TIME DEPENDENT 
PROFILE 

The form (|18|) suggests that beyond the instabil- 
ity the optimal profile is a fixed shape moving at a 
constant velocity v 

p = g(x — vt) . 

Due to conservation law © the current is then 

t) = jo - vp + vg(x - vt) 

If such a profile is the optimal profile, then the vari- 
ational principle © reduces to 



GUo) 



dx 



g{x),vj 2a(g(x)) 



jo - vp + vg{x) 



-D{g{x))g'{x)-vcj(g{x)) (19) 



This is of the form (the term linear in g' gives a 
null contribution due to the periodic boundary con- 
ditions) 



G(jo) = -inf / dx[X(g)+g' 2 Y(g)] (20) 



where 



and 



X(g) 



[jo - vpo + vg- ^cr(g)] 2 
D 2 (g) 



Y(g) 



The optimal v in l|19|l is then given by 

J dx (9 - po)U l , - U(,{9)) fdxt^Bil 
- = -jo- 



(21) 

this last simplification being due to the constraint 
J g(x)dx = p . With this constraint and for a fixed 
v, a variational calculation of the optimal g in l|20() 
shows that g should satisfy 

X\g) - 2Y{g)g" - g' 2 Y'{g) = C 2 

Multiplying both sides by g 1 allows one to integrate 
once so that g satisfies 



X(g) - g' 2 Y(g) =d + C 2 g 



(22) 



where C\ and C 2 are constants (which is an exten- 
sion of the equation (15) of 7] to the case of the 
ring). 

For fixed jo, po and i>, if one denotes by g\ and 
32 the two extrema of the profile g (generically, the 
profile g{x) is a periodic function of period 1 with 
a single minimum g\ and a single maximum g 2 ), 
one can determine the constants C\ and C 2 by 1221) 
in terms of g\ and 172 (as X(gi) = C\ + C 2 g\ and 
X(g 2 ) = C\ + C 2 g 2 ). The differential equation 
determines the whole profile (up to a translation on 
the ring) and the constants g\ and g 2 are then fixed 
by the fact that 



■ {si) 



dx 



92 dg_ 

91 9' 




Y{g) 



X{g) C 2 g 



dg 



and 



2 



fi-i 



Y(g) 



X{g) C 2 g 



dg 



IV. EXACT NUMERICS FOR THE 
WEAKLY ASYMMETRIC EXCLUSION 
PROCESS ON A RING 

We wrote a program to calculate exactly p(X) for 
the weakly asymmetric exclusion process (WASEP) 
on a ring of sites with P = N p particles. In the 
simple symmetric exclusion process (SSEP), each 
particle jumps to its right at rate h and to its left 
at rate \ and the functions D(p) and a(p) are given 
1 by 



-Dssep — 2 



(23) 
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WASEP: asymmetry v 
T 



10 and density 1/2 
r 



A^u(A) 




Figure 1: /i(A) as denned by 111) II for the weakly asym- 
metric exclusion process on a ring of N = 6, 10, 14, 18, 22 
sites when the asymmetry parameter v = 10 and the 
density 1/2 (thin lines). As N increases, the data in- 
crease and seem to accumulate to the value predicted 
(dashed line) by assuming that the optimal profile is the 
one discussed in section fTTTI If the optimal profile had 
been flat, the curve would have been predicted by 1131 
(thick line). The horizontal dotted line gives the value 
of \x below which the flat profile becomes unstable 1251 . 



CSSEP = p(l - p)- 



(24) 



If one introduces a weak electric field to the right, 
the model becomes the WASEP and the rates be- 



come 



, , to the right and | - -j^. 
As the evolution is a Markov process, one can 
build, as explained in |18lll9| from the Markov ma- 
trix, a A-dependent matrix, the largest eigenvalue 
of which is /i(A) defined by Q10[l. According to the 
linear stability analysis, the flat profile becomes un- 
stable {HO) for 

j 2 < p{l - p)[v 2 p{l - p) - it 2 } 

or by lfl7| for 



A^flat(A) < - — . 



(25) 



We have calculated the exact eigenvalue p(X) for 
lattice sizes from N = 6 to 22, at density 1/2 for 
an asymmetry v = 10 (in order to avoid negative 
rates for small system sizes, we have replaced in 
our programs the rates \ ± ^ by exp[±f/JV]/2). 
The results show a rather quick convergence with 
increasing N towards the value corresponding to a 
shape determined by H19|) . Clearly the flat profile 
gives a value too low, incompatible with the numer- 
ical data. 



V. A BRIDGE BETWEEN A WEAK AND A 
STRONG ASYMMETRY 

Assuming that the optimal profile is the one dis- 
cussed in section IIIII if one tries to make v large, 
and one writes 



Jo = viq 



(26) 



the moving profile which satisfies 11221) becomes very 
steep in the regions where it varies, and it takes the 
form of a step function with two constant values g\ 
and g 2 separated by two discontinuities 



g(x) =gi for < x < y 
g(x) = .92 for y < x < 1 



(27) 



so that the parameters g\ , g 2 and y are related to po 
and «o by 



Po = ygi + (1 - v)g2 



(28) 



«o = ya(gi) + (1 - y)(r(g 2 ) . (29) 
The expression of the velocity l|21|) then becomes 
c(ff2) - cr(gi) 



.92 - 5i 



(30) 



Then using the fact that g' vanishes when g(x) — 
,9i or g 2 and replacing (|26I28I29I30() into (J2H implies 
that the constants C\,C% (in (|22Jl ) vanish at order 
v 2 . Thus asymptotically in v, one can rewrite 1221) 
as 



1 



9' 



D(g) 2 ( 9l 



92) 



32KS1) -<<?)) (31) 



+5iO(5) ~ cr(g 2 )) + g(cr(g 2 ) - a(gi)) 



and 



G(jo) 



dx 



o{g){gi - S2) 2 



g2{a{gi) - <j(g)) 



+gi(a(g) - a{g 2 )) + g(a(g 2 ) - a( gi ))\ . (32) 

If one replaces g by its expression (|2"7|l , one gets that 
the order v 2 of G(jo) vanishes. The next order in the 
large v expansion is dominated by the rounding-off 
of the discontinuities in l|27(l as given by 1)31(1 . As the 
profile g(x) is composed of two monotonic parts one 
can then use H31|l into (|32|) and obtain for g 2 > 91 



D(g) 



G(io) = -2v / dg- 

J gi (.92 - .9iM.9) 



c r (5)(32 - .9i) 



-c(5i)(52 ~ .9) ~ cr(f2)(.9 - .9i) 



G 



It is remarkable that y is not present in this expres- 
sion. In the case of the weakly asymmetric exclusion 
process on a ring, expressions (|23I24|) of D(p) and 
a(p) lead for large v to 



G(jo) 



92 -9i -9\92 In 



92 

9i 



(33) 



-(l-5i)(l-S2)ln 



l-9i 



1 - 92 

In this case (|3()l) becomes u = — fli — 32)- 



If one takes formally v — N, the hopping rates 
1/2 ± v/2N become 1 and 0, so that the model re- 
duces to the totally asymmetric exclusion process 
and one gets from H7I26I33[) 



Pro 



Or 
T 



la 



exp 



-T 



92 



-(l-ffO(l-fla)ln 



91 - gm in 

1-51 



92 



1 



92 



91 
(34) 



As we will see it in section fvTl this is exactly the 
large deviation function predicted by the Jensen- 
Varadhan theory to maintain a profile l|27(l 
formed of a shock and an antishock in the totally 
asymmetric exclusion process. Other aspects of the 
relation between the large deviation functional of 
the weakly asymmetric exclusion process and the 
Jensen- Varadhan functional in systems with open 
boundary conditions will be presented in |2l|. 



VI. LARGE DEVIATIONS OF THE 
CURRENT IN THE TOTALLY 
ASYMMETRIC EXCLUSION PROCESS 

In the totally asymmetric process, each particle 
jumps to its neighboring site, on its right, at rate 
1, if the target site is empty (and there is no other 
jump). 

The large deviation function of the current of the 
totally asymmetric exclusion process on a ring of N 
sites, with P particles, has been calculated exactly 
[l8L |2^ . If Qt is the total number of jumps during 
time T over a given bond on the ring, one knows 
that for large T, 



(e XQT ) 



(35) 



and explicit expressions of /i(A) has been obtained 
for all N and P by the Bethe ansatz [T l l ^ . 

For large TV, it was shown in particular (equa- 
tion (53) of [22] with the proper redefinition of the 
parameters) that for A < 



MA) = - 



(1 



)(1 



iHi-po) 



(36) 



We are going now to argue that this result can be 
understood, by assuming that l|35|) is dominated by 
configurations of the form (|27l) moving at a velocity 
v = 1— fli— fl2- These density profiles are everywhere 
constant except for a shock (at some position z with 
g(z — 0) = g\ and g(z + 0) = 32 for 32 > fli) and an 
antishock (at position z + y with g(z + y — 0) = 32 
and g[z + y + 0) = 31). From the Jensen- Varadhan 
theory |20|. the probability of maintaining such a 
shape moving at this velocity on the ring over a very 
long period of time T reduces to the probability of 
maintaining an antishock between the densities 32 
and <?i moving at velocity v. The probability of the 
latter event, which we denote by Pt(<?i, 32) is given 
E3 by 



-T 



-Pt(3i,32) ~ exp 

-(l- gi )(l- 92 )\n 



92 - 9i - 9W2 In 



92 

9i 



l-9i 



(37) 



,1-32 

The corresponding integrated current Qt is 
Qt = T[y 9l (l - 9l ) + (1 - 3)32(1 - 32)] 



since over a long period of time T a given bond 
spends a fraction y of the time at density 31 and 
1 — y at density 32 . 

Therefore, if the configurations of the form l|27(l 
dominate the large deviations of the current, one 
expects 



//(A) = max {A[y3i(l - 31) + (1 

y,gi,g2 



3)32(1 



32 "3i -3i32 In |f- - (1 - fli) (1 -P2) In 



52)]- 

(38) 

where the maximum has to satisfy the constraint 
Po = V9i + (1 - 3)32 (39) 

A calculation of the optimum in 138|) . with the 
constraint l|39l) leads to 



o A(i-p ) 



9i 



1 



32 



and (|3*5|l becomes l{3!j|) . 

This shows that the result of the Bethe ansatz 
(|36|l can be physically understood in terms of an op- 
timal profile which takes the form of the step func- 
tion (|27|l . The probability of maintaining this pro- 
file is given by the Jensen- Varadhan expression l|37(l 
which in fact is identical to (|34|l obtained, in the 
large v limit, for the WASEP from the hydrody- 
namic large deviation theory. 

That the fluctuations are due, in the strong asym- 
metric case, to configurations formed by a gas of 
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shocks and antishocks has been already pointed out 
by Fogedby [H HI HJ . The calculation of this sec- 
tion shows that the large deviation of the current, in 
the range A < 0, can be understood quantitatively in 
terms of a single pair of shock-antishock. Whether 
the Fogedby theory would allow to understand all 
the current fluctuations, including the range A > 
where the expression of //(A) is more complicated 
|22| than i(36() . remains an interesting open question. 



VII. CONCLUSION 

In the present work we have determined the limit 
of stability 1(16117(1 of a flat profile for a diffusive 
lattice gas on a ring. This instability beyond which 
the optimal profile becomes modulated is of the 
same nature as the phase transition found for several 
other non-equilibrium systems [SIISEIIIII- As the 
calculation is based on a local stability analysis, one 
cannot of course exclude first order transitions, i.e. 
that the flat profile might become globally unstable. 

In section IIVI we have obtained numerical evi- 
dence that the macroscopic fluctuation theory pre- 
dicts correctly the large deviation function of the 
current for the weakly asymmetric exclusion pro- 
cess. The numerical results are consistent with the 
second order phase transition predicted in section 
ITT1 and with a modulated density profile moving at a 
constant velocity as suggested in section ITTT1 These 
results could in principle be confirmed by solving 
the Bethe ansatz equations for the WASEP, since 
ju(A) can be calculated exactly for the ring geome- 
try mill. 

It would be interesting to extend the results of 
the present work to the case of open boundary con- 
ditions. One difficulty is that the time-independent 
profile, found in 0, is much more complicated than 
the flat profile for the ring geometry, and we did not 
succeed so far to obtain the condition which would 
generalize 1(16117(1 for this open geometry. 

Lastly, we noticed that the large deviation func- 
tion ((34(1 obtained for the weakly asymmetric dif- 
fusion process in the large drift limit is identical 
to the one predicted for a strong asymmetry by the 
Jensen- Varadhan theory (see also HH). 

Despite this bridge between the large deviation 
function of the current of weakly and strongly asym- 
metric systems, and some recent results on zero- 
range processes 0, a theory of current fluctua- 
tions for strongly asymmetric lattice gas such as the 
ASEP with open boundary conditions remains an 
open problem. 



VIII. APPENDIX: A DERIVATION OF lHl5l) 

We present here an heuristic derivation of the hy- 
drodynamic large deviations Let us consider a 
system of N sites and decompose it into N/l boxes 
of / sites each. Let us define the density pi (t) in box 
i at time t and qi(t) the total number of particles 
transfered from box i to box i + 1 during a time in- 
terval t,t + r (this time r should be large enough for 
the qi{t) to be a Gaussian characterized by its av- 
erage and its variance as in (lll-il) . but short enough 
compared to the characteristic time of variation of 
the densities Pi(t)). 

If one writes that the qi(t) are Gaussian, one gets 



Pm(qi(s),pi(s)) ~ exp 



T/t N/l 



EE 2.( P !( S ))r (g'( g ) + 



k=l i=l 



where k — sr. The factor v /N comes from the weak 
asymmetry of the jumps. Clearly the conservation 
of the number of particles gives 

p t (s + T)= pi (s) + (40) 



Now if one takes a continuous limit by writing 

I s \ 



Pi(s)=p{i N , N2 j 
and one defines a rescaled current by 



(41) 



one gets 

Pro(j(x,s),p(x,s)) ~ exp 

[j(x, 8) + D(p(x, 8))&g£ Mp(x, s)) 



(42) 



-iV" 1 / ds I dx 
o Jo 

2 



2a(p(x,s)) 



which is exactly (1 lloll . Furthermore (|40|) with the 
scaling 141142(1 leads to ©• 

Acknowledgments: We thank C. Bahadoran, L. 
Bertini, A. De Sole, D. Gabrielli, G. Jona-Lasinio, 
C. Landim for very helpful discussions. 



8 



[1] L. Bertini, A. De Sole, D. Gabrielli, G. Jona- 
Lasinio, C. Landim, Fluctuations in stationary non 
equilibrium states of irreversible processes, Phys. 
Rev. Lett. 87, 040601 (2001) 

[2] L. Bertini, A. De Sole, D. Gabrielli, G. Jona- 
Lasinio, C. Landim, Macroscopic fluctuation the- 
ory for stationary non equilibrium states, J. Stat. 
Phys. 107, 635-675 (2002) 

[3] C. Kipnis, S. Olla, S. R. S. Varadhan, Hydrody- 
namics and large deviations for simple exclusion 
processes, Comrnun. Pure Appl. Math. 42, 115-137 
(1989) 

[4] H. Spohn, Large scale dynamics of interacting par- 
ticles, Springer (1991) 

[5] C. Kipnis, C. Landim, Scaling limits of interacting 
particle systems, Springer (1999) 

[6] L. Bertini, A. De Sole, D. Gabrielli, G. Jona- 
Lasinio, C. Landim, Current fluctuations in 
stochastic lattice gases Phys. Rev. Lett. 94, 030601 
(2005) 

[7] T. Bodineau, B. Derrida, Current fluctuations in 
nonequilibrium diffusive systems: An additivity 
principle Phys. Rev. Lett. 92, 180601 (2004) 

[8] L. Bertini, D. Gabrielli, J. L. Lebowitz, Large 
deviations f or a s tochastic model of heat flow, 
|cond-mat/050168l1 

[9] C. Kipnis, C. Marchioro, E. Presutti, Heat-flow in 
an exactly solvable model J. Stat. Phys. 27, 65-74 
(1982) 

[10] B. Derrida, J. L. Lebowitz, E. R. Speer, Free energy 
functional for nonequilibrium systems: an exactly 
solvable case, Phys. Rev. Lett. 87, 150601(2001). 

[11] B. Derrida, J. L. Lebowitz, E. R. Speer, Large de- 
viation of the density profile in the symmetric sim- 
ple exclusion process, J. Stat. Phys. 107, 599-634 
(2002) 

[12] C. Enaud, B. Derrida, Large deviation functional of 
the weakly asymmetric exclusion process, J. Stat. 
Phys. 114, 537-562 (2004) 

[13] B. Derrida, B. Dougot, P.-E. Roche, Current fluctu- 
ations in the one-dimensional symmetric exclusion 
process with open boundaries, J. Stat. Phys. 115 
717-748 (2004) 

[14] F. van Wijland, Z. Racz, Large deviations in weakly 
interacting boundary driven lattice gases, J. Stat. 
Phys. 118, 27-54 (2005) 

[15] R. J. Harris, A. Rakos, G. M. Schiitz, Current fluc- 
tuations in the zero-range process with open bound- 
aries, cond- mat/0504711 

[16] S. Pilgram, AN. Jordan, E.V. Sukhorukov, M. 
Buttiker, Stochastic path integral formulation of 
full counting statistics, Phys. Rev. Lett. 90, 206801 



(2003) 

[17] A.N. Jordan, E.V. Sukhorukov, S. Pilgram, Fluctu- 
ation statistics in networks: A stochastic path inte- 
gral approach, J. Math. Phys. 45 4386-4417 (2004) 

[18] B. Derrida, J.L. Lebowitz, Exact large deviation 
function in the asymmetric exclusion process, Phys. 
Rev. Lett. 80, 209-213 (1998) 

[19] P.-E. Roche, B. Derrida, B. Dougot, Mesoscopic full 
counting statistics and exclusion models, Eur. Phys. 
J. B 43, 529-541 (2005) 

[20] L. Jensen, The asymmetric exclusion process in 
one dimension, Ph.D. dissertation, New York Univ., 
New York, 2000 

S.R.S. Varadhan, Large deviations for the asym- 
metric simple exclusion process, Stochastic analysis 
on large scale interacting systems, 1-27, Adv. Stud. 
Pure Math., 39, Math. Soc. Japan, Tokyo, (2004) 

[21] T. Bodineau, B. Derrida, Current large deviations 
for asymmetric exclusion processes, in preparation 

[22] B. Derrida, C. Appert, Universal large deviation 
function of the Kardar-Parisi-Zhang equation in one 
dimension J. Stat. Phys. 94, 1-30 (1999) 

[23] I. Vilfan, R. K. P. Zia, B. Schmittmann, Sponta- 
neous structure formation in driven systems with 
2 species - exact-solutions in a mean-field theory 
Phys. Rev. Lett.73, 2071-2074 (1994) 

[24] G. Fayolle, C. Furtlehner, Dynamical windings of 
random walks and exclusion models. Part I: Ther- 
modynamic limit in Z(2) J. Stat. Phys. 114, 229- 
260 (2004) 

[25] M. Clincy, B. Derrida, M.R. Evans, Phase transi- 
tion in the ABC model Phys. Rev. E 67, 066115 
(2003) 

[26] H.C. Fogedby, Morphology and scaling in the noisy 
burgers equation: Soliton approach to the strong 
coupling fixed point Phys. Rev. Lett. 80, 1126-1129 
(1998) 

[27] H.C. Fogedby, Soliton approach to the noisy burgers 

equation: Steepest descent method Phys. Rev. E 

57, 4943-4968 (1998) 
[28] H.C. Fogedby, A. Brandenburg, Solitons in the 

noisy Burgers equation Phys. Rev. E 66, 016604 

(2002) 

[29] D.S. Lee, D. Kim, Large deviation function of the 
partially asymmetric exclusion process Phys. Rev. 
E 59, 6476-6482 (1999) 

[30] D. Kim, Bethe-ansatz solution for crossover scal- 
ing functions of the asymmetric XXZ chain and the 
Kardar-Parisi-Zhang- type growth-model Phys. Rev. 
E 52, 3512-3524 (1995) 



